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<D Abstract 

Dynamical chiral symmetry breaking (D^SB) is studied in the Nambu- 
Jona-Lasinio model for an arbitrary combination of external constant electric 
and magnetic fields. In 3 + 1 dimensions it is shown that the critical coupling 
constant increases with increasing of the value of the second invariant of elec- 
tromagentic field E ■ B, i.e. the second invariant inhibits D\SB. The case 
of 2 + 1 dimensions is simpler because there is only one Lorentz invariant of 
electromagnetic field and any combination of constant fields can be reduced 
to cases either purely magnetic or purely electric field. 
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In works [1-3] it was first shown in the so called ladder approximation that 
quantum electrodynamics (QED) in the regime of strong coupling has a new phase 
with dynamically broken chiral symmetry. The new phase of QED possesses very 
interesting properties from the theoretical viewpoint [4-6]. There were attempts to 
use them for modelling electroweak symmetry breaking in technicolor-like models 
[7]. However, at present the new phase of QED has little relevance to experiment 
because we need a strong coupling constant a c ~ 1 (recall that the physical value 
of electromagnetic coupling constant is a c = <C 1). 

However, as was suggested in [8, 9], the situation may drastically change in the 
presence of strong external electromagnetic fields where dynamical chiral symmetry 
breaking (DxSB) may occur at the regime of weak coupling. A breakthrough in this 
direction was made in [10, 11], where, in the framework of the Nambu-Jona-Lasinio 
(NJL) model [12], it was shown that DxSB takes place in an external constant 
magnetic field at any small attraction between fermions both in 2 + 1 and 3 + 1 
dimensions (note that the fact that external magnetic field enhances DxSB was 
first noted in the NJL model in [13] (see also [14])). 

As shown in [10, 11], in the infrared, the dynamics of fermions in magnetic field 
in 2 + 1 and 3 + 1 dimensions resembles the dynamics of fermions in + 1 and 1 + 1 
dimensions, respectively. Therefore, we have an effective reduction of dimension of 
space-time by 2 units and as a result the critical value of the coupling coustant in 
external magnetic field is equal to zero. It was latter shown in [15] that the same 
effect takes place in QED in external magnetic field. Note that although the critical 
value of coupling constant is zero extremely strong magnetic fields (\B\ > 10 13 G) 
are necessary for experimentally significant consequences because the correction to 
the physical mass of electron is very tiny for weak magnetic fields. 

The case of constant electric field was considered in [13] where it was shown 
that the value of the critical coupling constant is more in this case than in the case 
without electric field. In the present work we study D%SB in the NJL model in the 
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case of an arbitrary combination of constant electric and magnetic fields in 3 + 1 
and 2 + 1 dimensions. We first consider the case of 3 + 1 dimensions. As well 



known, electromagnetic field has two Lorentz invariants f\ = \F ilv F iil ' = B 2 — E 2 
and / 2 = \s^ val3 F^ v F a f3 = E • B. Since D% SB was already studied in cses of purely 
electric and magnetic constant external fields when only the first invariant f\ of 
electromagnetic field is not equal to zero, in the present work we consider D%SB in 
the case where f 2 ^ 0. 

The Lagrangian of the NJL model [12] in an external electromagnetic field reads 

N s~i N 

3=1 Z 3=1 

where is the covariant derivative = <9 M + ieA^ and j = 1, 2, . . . , N flavor 
index. Lagrangian (1) is invariant with respect to the U L (N) x U R (N) chiral group. 
By using auxiliary fields tt and a, we can rewrite (1) in the following form: 

N r 

tf Vn nr.. xb irv i ;- -~.nl/. _ _L (a 2 + 7lj 



r - l 

£ = E Mj-fD^j - M>> ; + /-r,;T / )vl/ ; - — (a 2 + tt 2 



(2) 



By taking integrals over fermion fields, we obtain the effective action for tt and 
a fields 

jv , 
tt) = -< £ Tr Ln [i^D, - (a 3 + ^vr,)] -—J d 4 x(a 2 + ttJ). (3) 

To obtain the effective potential for a and tt fields, it suffices to consider the case 
of constant fields a = const, tt = const. Since the effective action is invariant with 
respect to the Ul(N) x Ur(N) chiral symmetry, the effective potential depends on 
tt and a fields only through the chirally invariant combination p 2 = X^LiO^f + )• 
Therefore, in what follows it is sufficient to set tt^ = , o~k — for k = 2, . . . ,N 
and consider the effective potential only for the field U\ which we simply denote a. 
Thus, 

r(<j) = -*Tr Ln [i^D, - a] - ± J d 4 xa 2 . (4) 
By using the method of proper time [16, 17], we represent the first term in (4) as 
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follows: 

- iTrlmiiDrf - a) = — % -TrLn[D 2 + a 2 ) = [ ^-tr(x\e- is{D ' 2+a2) \x)dsd 4 x (5) 

As well known [17], vacuum of QED is not stable in an external electric field 
and the effective potential has an imaginary part which defines the rate of birth 
of fermion-antifermion pairs from vacuum per unit volume. Since we study the 
problem of DxSB, we can ignore this effect and consider only the real part of 
effective potential which is equal to 

a 2 N r°° 1 2 

Via) = -r 7 z + Tr^v.p. / ds-e~ sa M coth(Ms)L cot(Ls), (6) 
2G an 2 Ji/a 2 s 

where L 2 = e 2 ^f~ f \ M 2 = e 2 ^f +} \ and f 1 = B 2 - E 2 and f 2 = E • B 
are two invariants of electromagnetic field. In (6) we introduced a cut-off and 
v.p. of the integral in s is present because we consider only the real part of the 
effective potential (recall that the imaginary part of the effective potential is given 
by residues in poles of cot(Ls)). The gap equation SV/Sa\ a=m = has the form 

1 N 

G An 2 "' 1 " A/A2 

This gap equation was investigated in cases where only the first invariant of 
electromagnetic field is not equal to zero, i.e. for cases of purely electric and magnetic 
external fields. In this paper we study how the presence of nonzero electric field 
parallel to magnetic field (E ■ B ^ 0) affects DxSB. By using some inequalities, we 
first analytically obtain an estimate from below for the critical coupling constant. 
We add and subtract 1/s to Mcoth(Ms) in the gap equation (7). Then 

v.p. / dse~ sm Mcoth(Ms)Lcot(Ls) = A 2 L + 

Jl/A 2 2 

r°° 2 

v.p. j dse~ sm (Mcoth(Ms) - l/s)L cot(Ls), (8) 

where we used the result [13] v.p. J^ A2 dse~ sm2 L^±^ = A 2 - f L. Further, we 
represent the integral in (8) as a sum of two integrals Jy A 2 = j^ L A 2 + /J? (note that 



r°° 2 

i.p. ^ ^ dse~ sm M coth(Ms)L cot(Ls) = 0. (7) 



is the first zero of cot(Ls)). We now consider the integral from ^ to infinity. 
Since cothx < 1/x + 1 for x > 0, we have 

P OO P oo 

v.p. I dse- sm2 (M coth{M s) -1/ s)L cot(Ls) <v. p. I dse- sm2 ML cot(Ls). (9) 
Integrating by part, we obtain 

f°° 2 f°° 2 

v.p. I dse~ sm MLcot(Ls) = M I e~ m s rf(ln(2|sinLs|)) 

= -Mm(2|sin-|) + m 2 M / e~ m Tn(2|sinLs|)ds. (10) 

By using the formula [18] 

7T 

and the fact that the integral M ln(2|sinLs|)ds is finite, we conclude that the 
integral m 2 M e~ m Tn(2|sinLs|)ds tends to zero on the critical line (where m 2 — > 
0). Thus, v.p. dse- sm2 MLcot{Ls) = -Mln2. Therefore, in view of (9), we have 



f°° 2 

v.p. I dse~ sm (Mcoth(Ms) - l/s)Lcot(Ls) < -Mln2. (12) 
•'ye 

TV 

It remains to estimate from below the integral j^h 2 <is(Mcoth(Ms) — l/s)L cot(Ls). 
We have jf A2 ds(Mcoth(Ms) - l/s)L cot(Ls) < (M cot (Ms) - l/s)f (because 
cotx < 1/x for x in the interval from to |). If M 3> L, then we use the estimate 
cothx < 1/x + 1 because coth(Ms) is approximately 1 near the upper limit of 
integration. Therefore, in this case 

[ 2L ds(Mcoth(Ms) - lfs)Lcot(Ls) < Mln^— . (13) 
A/a 2 2L 

If M <C L, then coth(Ms) <C 1 in the interval of integration. By using the estimate 
cothx < i + x/3 (| and x/3 are simply two first terms of the Teylor expansion of 
cothx), we obtain 

[ 2L ds(Mcoth(Ms) - l/s)Lcot(Ls) < -z-^ - 1/A 2 ). (14) 
Ji/A 2 3 2L 



We now analyse the obtained results. We assume in what follows that |/i| 3> I/2I. 
In the magnetic-type case (j\ > 0), we have L M^) 1 / 2 and M m \e\f^ 2 . By 
using (8), (12), and (13), we obtain the following estimate from below for the critical 
coupling constant in the magnetic-type case: 

1 - 27^ + A^ ln 2l - A^ lni l + |e|^ln^7^ 

where g cr is the dimensionless critical coupling constant g cr = 47r2 ^ A2 . It directly 
follows from (15) that the presence of electric field parallel to magnetic field is very 
important. Indeed, if / 2 7^ 0, then g cr is no longer equal to zero (even if the magnetic 
field is very strong \B\ ~ A 2 ) in contrast to the case of purely magnetic field where 
g cr = 0. If f 2 increases, g cr is also increases. If f 2 — > 0, then g cr — > 0, i.e. we recover 
the result obtained by Gusynin, Miransky, and Shovkovy [10]. In the electric-type 
case (fx < 0) we have L « Hl/il 1 / 2 and M « l e l(|^| ) 1/2 - % usin S ( 8 )> ( 12 )> and 
(14), we obtain 

9 " - 1 - S + iW- 1/A 2 ) - K 1 - *|p : g(^jv^- (16) 

It follows from (16) that in this case g cr for f 2 7^ is more than g cr = * 1/2 [13] 

1 ^ I e 1 1 J 1 1 ' 

1 iX^ — 

in the case f 2 — 0. As / 2 goes to zero, our estimate coincides with the result obtained 
by Klevansky and Lemmer [13] in the electric-type case. Thus, we conclude from 
the obtained estimates that the second invariant of electromagnetic field inhibits 
DxSB. In magnetic-type case it looks rather natural (indeed, if f 2 7^ 0, then it 
means that E 7^ and we know that electric field inhibits DxSB). However, in 

— * — * 

the electric-type case it appears unlikely. Indeed, let first £^0,5 = (therefore, 
f 2 — 0). It is natural to assume that g cr should decrease in the case f 2 7^ because if 
f 2 7^ it means that B 7^ and we know that magnetic field assists DxSB. We can 

— * — * 

understand the cause of growth of g cr with increasing of E ■ B as follows. Since we 
study the dependence on the second invariant, we keep the first invariant B 2 — E 2 
unchanged. Without loss of generality we can assume that E || B (if not, one can 
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perform an appropriate Lorentz transformation). If we increase f 2 , then in order to 

— * — * 

keep the first invariant unchanged we have to increase both B and E. Therefore, 
there is a competition between increasing of B and increasing of E. It turned out 
that qualitatively increasing of E is more significant for g cr than increasing of B for 
any fi, therefore, g cr always grows with increasing of f 2 . 

We found a rather rough analytic estimate from below for the critical coupling 
constant. To obtain a more accurate dependence of the critical coupling constant on 
f 2 , we numerically calculate the integral in (7). A typical dependence of g cr on f 2 
in the electric- type case is shown in Fig.l (this figure corresponds to = — 10~ 4 ) 
and in the magnetic-type case in Fig.2 (where -p = 1CT 4 ). 

We see from these figures that the critical coupling constant increases with in- 
creasing of f 2 for any value of f\. In the electric- type case g cr is always more that 
1. In the magnetic-type case g cr abruptly drops to zero as f 2 tends to zero. We 
also numerically calculated g cr in the case where the first invariant is zero f\ = 
(\E\ = \B\) and obtained a dependence which is similar to the electric-type case, i.e. 
g cr increases with increasing of f 2 . Thus, the numerical analysis of the gap equation 
confirms that the second invariant of electromagnetic field inhibits D\SB. 

We now consider the case of 2+1 dimensions. We use the reducible 4-dimensional 
representation of the Dirac algebra for fermion field in order that the model possess 
a chiral symmetry (in fact, there is two chiral symmetries with 75 and 73 matrices, 
for more details see [19]) and we do not study parity breaking. Thus, we have the 
following gap equation for parity conserving mass: 

m = 2iGtr{S<x l) {x,x)), (17) 

where is the fermion propagator in an external constant electromagnetic field 
which has the following form in the Fock-Schwinger proper time formalism: 

S { ™\x,x') = (id^-eA^-m)(-i) f drU A (x,x';r), (18) 

J— 00 
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where e 



and 

U A (x,x'; t) 



< x\ 



(id - eA) 2 - \o^F i 



x > 



m" 



87r 3 / 2 |r| 3 / 2 



exp 



re 



diA(i) + 



-{x- x')eF coth(eFr)(x - x') - 
-tr[\n ( )] + -<rFr + 

By taking trace over the Dirac indices in (13), we obtain 



m 



mG 



-itt/A 



TV 



o dr 



oo r 



3/2 



(eXr) cot(eXr), 



(19) 



(20) 



where X = y^F^F^ = y B 2 — E 2 (note that magnetic field is a pseudoscalar in 
2 + 1 dimensions, not axial- vector) . 

As well known, in 2 + 1 dimensions there is only one Lorentz invariant \F iW F ilv 
of electromagnetic field. Indeed, we can see from (20) that the gap equation de- 
pends on electromagnetic fields only through this Lorentz invariant. Consequently, 
an arbitrary combination of constant fields can be reduced to cases either purely 
magnetic or purely electric field. The case of constant magnetic field in the Nambu- 
Jona-Lasinio model in 2 + 1 dimensions was considered in [10]. Therefore, we study 
here only the case of constant electric field. 

In the case of external electric field the gap equation takes the form 



m 



mG 
mG 

7T 



-i7r/4 



—in/4 



/: 



00 ds 



,3/2 



e - sm + 



dr 



3/2 



r u c 
J -00 It 



3 ir(m 2 — ie) 



3/2 



[(eEr) coth(eFr) - 1], 



(21) 



where we explicitly wrote down the term which corresponds to the gap equation 
without external electric field and E — \E\. Further, by using the fact that 
rcothr — 1 = r(l + e2T 2 _ 1 — ^) and performing the change of variable x — |, for 



Rc 



i7r/4 jO dr r iT(m 2 -ie) 



|7p72 



(eEr) coth(eFr)] we get 



Re 



Im 



o) 



(22) 
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where Ie(i^%, \) = Jq° x a ~ 1 e l ~z~ x (^—r^j — ^j, [i\ = ^§ and we also used the equality 

ir^( cos{ T x) - shi( T x) ) =o - < 23) 

Ie^ei a ) ls an analytic function of a in the region < Re(a) < 2. 



For 1 < Re(a) < 2, by representing Ie{He) as a sum °f two integrals, we have 



7 B (/4, a) = -(^y-'ria - l)e^ + r(a)((a, I -A. (24) 



Performing an analytic continuation of (24) to the region < Re(a) < 2, we get 



2' v ^ /r ,v 2 
Consequently, the gap equation takes the form 

r3/ 2 .fflpfil 1 ^ 



tt 3 / 2 _ i 0F(e£) 



2GA 



2A 



Re 



— Im 



4 x -4> 



(25) 



+ <0 <*» 



From (26), we obtain the following value for the critical coupling constant: 

vr 3 / 2 1 g cr (0) 



9cr(E) 



2A l-o'Ca)® 1 / 2 l-a 2 (f; 



1/2 : 



(27) 



where g cr = °^0r is the dimensionless coupling constant in 2+1 dimension and 
a 2 = -f C(|) « 1-29. 

Thus, the presence of external constant electric field increases the value of critical 
coupling constant. Note that the same is true in 3 + 1 dimensions (see [13]). It is 
not difficult to show that in the vicinity of critical point 



m 



C (E cr — E) 



- RW2 



(28) 



where C 2 = -2^| ~ 1.12. 

Thus, this phase transition is a phase transition of the second order. 

Since we have only one Lorentz invariant of electromagnetic field X 2 = -^F^F^ = 
B 2 — E 2 in 2 + 1 dimensions, by using an appropriate frame, the general case of 
non-zero constant electromagnetic field can be reduced to the cases of purely electric 
or purely magnetic fields. 
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Figure captions 

Figure 1: The dependence of the critical coupling constant on the second invariant 
in the magnetic-type case. 

Figure 2: The dependence of the critical coupling constant on the second invariant 
in the electric-type case. 
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